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QUANTUM  CORRECTIONS  TO  THE  THERMODYNAMIC  PROPERTIES 
OF  LIQUIDS,  WITH  APPLICATION  TO  NEON 


By  O.  K.  Rice* 


ABSTRACT 

The  quantum  correction  has  been  expressed  as  a  deviation  rj  from  RT  In  Q^,  where  Q^  is  the  clas¬ 
sical  partition  function.  Changes  in  volume,  energy,  specific  heat,  etc.,  caused  by  quantization  can  be 
expressed  in  terms  of  r}  and  certain  thermodynamic  quantities,  r?  itself  is  expressed  in  terms  of  ther¬ 
modynamic  quantities.  The  equations  obtained  give  for  solid  neon  the  same  results,  within  the  limits 
of  error,  as  the  usual  treatment  of  a  Debye  solid  with  the  aid  of  the  Griineisen  equation.  To  apply  the 
equation  to  liquid  neon,  various  thermodynamic  quantities  have  been  obtained  from  known  values  for  ar¬ 
gon  by  use  of  the  law  of  corresponding  states.  The  difference  in  volume  and  specific  heat  of  liquid  neon 
from  the  values  expected  from  the  law  of  corresponding  states  has  been  computed  and  found  to  agree 
reasonably  well  with  the  observations.  However,  the  melting  point  and  the  AE  and  AS  of  fusion  are  found 
experimentally  to  be  much  closer  to  the  classical  values  than  one  would  expect  from  the  theory.  No 
explanation  of  this  discrepancy  is  apparent. 


+  *  *  ♦  + 


It  has  frequently  been  noted  that  the  properties  of  neon  are  somewhat  out  of  line  with  those  of  the 
other  condensed  rare  gases.  These  deviations  have  been  ascribed  to  quantum  effects  in  neon,^  which, 
it  is  believed,  would  otherwise  obey  the  law  of  corresponding  states.  Quantum  effects  in  the  solid  are 
easily  taken  care  of  with  the  help  of  the  Debye  theory;  the  liquid  presents  more  of  a  problem. 

The  quantum  corrections  for  a  liquid  have  recently  been  given  in  a  general  form  by  Mayer  and 
Band.^  Their  results,  however,  involve  some  multiple  integrals  containing  distribution  functions  and 
second  derivatives  of  the  interatomic  potential  function.  Therefore,  they  are  difficult  to  evaluate,  and 
it  seems  desirable  to  derive  some  rougher  equations  which  will  be  more  easily  handled.  It  will  be  the 
purpose  of  this  paper  to  express  the  magnitude  of  the  quantum  effects  in  terms  of  thermodynamic  quan¬ 
tities  which  have  either  been  measured  or  may  be  approximated  from  the  law  of  corresponding  states, 
and  apply  the  results  to  neon. 

1.  EQUATION  FOR  A  ONE -PHASE  ISOTHERMAL  AND  ISOBARIC  SYSTEM 

We  shall  express  the  quantum  deviation  by  means  of  a  quantity  17  defined  by  the  equation t 

In  Q  =  In  Q^  +  T] 

Here  Q  is  the  actual  partition  function,  and  Qc  is  the  classical  partition  function  of  the  system,  at  the 
given  temperature  T  and  molal  volume  V.  Q^  is  found  by  taking  the  limiting  form  of  Q  for  high  temper¬ 
atures,  and  inserting  the  actual  value  of  T  in  this  limiting  form.  In  calculating  Q,  only  the  thermal 
part  of  the  energy  E^  and  the  zero -point  energy  E2;  are  considered.  The  Helmholtz  free  energy  is  there¬ 
fore  given  by 

A=“RTlnQ+Ep  (2) 


*On  leave  1946-47  from  the  University  of  North  Carolina,  Chapel  Hill,  North  Carolina, 
tif  n  is  very  small  we  may  write  equation  1  in  the  form  Q  =  Qc(l  +77),  which  brings  out  the  signi¬ 
ficance  of  Tj.  However,  to  actually  define  rj  by  this  equation  would  be  less  convenient  for  purposes  of 
computation  than  to  use  equation  1. 
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where  Ep  is  the  potential  energy  (excluding  any  potential  energy  involved  in  Ej.  and  E^);  that  is,  Ep  is 
the  value  the  energy  would  have  at  absolute  zero  with  all  the  atoms  at  their  positions  of  equilibrium. 
For  a  solid,  the  positions  of  equilibrium  form  a  regular  lattice,  but  this  is  of  course  not  true  for  a 
liquid.  From  equation  1  this  may  be  written 

A  =  -RT  In  Qc  -  RT  T]4-  Ep  (3) 

=  -  RTt] 

In  this  equation  we  assume  either  that  Ep  is  the  same  for  the  classical  and  quantum  case  if  the  volume 
remains  fixed,  or  better,  that  3.ni7  change  is  included  in  r].  Since  Ep  is  the  energy  after  subtracting  out 
all  thermal  and  zero -point  energy,  it  would  appear  that  any  change  in  it  which  could  be  caused  by  quan¬ 
tization  could  be  referred  to  a  change  in  the  average  coordination  number  caused  by  quantization.  Such 
a  change  could  occur  in  the  liquid,  but  may  be  included  in  r].  This,  then,  eventually  becomes  merged 
with  the  assumptions  regarding  t]  which  are  later  made,  in  order  to  estimate  it  for  the  liquid. 

The  Gibbs  free  energy  F  is  equal  to  the  Helmholtz  free  energy  if  the  pressure  p  =  0,  as  is  essen¬ 
tially  the  case  for  condensed  systems  around  atmospheric  pressure.  It  is  to  be  noted,  however,  that 
in  equation  3  A^  is  to  be  evaluated  at  the  same  molal  volume  as  A,  not  the  same  pressure.  If  A  is  eval¬ 
uated  at  zero  pressure,  must  be  evaluated  at  that  pressure,  5p,  at  which  the  system  would  have  the 
same  volume  if  the  classical  statistics  held. 

Since  the  pressure  is  given  by 

P  =  -{aA/av).p  (4) 

and  since  6p  is  the  pressure  the  system  would  have  if  w^ere  the  correct  partition  function,  we  may 
write 

Sp  =  -(aA(,/av)x  (5) 

The  true  pressure  being  zero,  we  obtain  from  equations  3,  4,  and  5 

6p  =  -RT(aJ5/8y)T  (S) 

If  we  let  be  the  classical  value  of  A  at  zero  pressure,  we  have* 

6p  '^P 

Ac^AcjO"^  J  ^^^c/'^P)q'  J 

G  ’  0 

Using  equation  5,  and  integrating,  assuming  (av/ap)rp  ^  to  be  constant, 

=  Ac,0  “  (5p/2)^  (0V/0p),j^p 

Aside  from  the  assumption  that  (av/ap)^  ^  is  constant,  this  equation  is  exact.  If  77  and  hence  6p  is  a 
small  quantity,  it  shows  that  A«  "  A^  q  is  a  small  quantity  of  the  second  order.  Such  quantities  may 
in  general  be  neglected  in  the  case  of  the  liquid,  but  may  not  be  entirely  negligible  for  the  solid.  The 
form  of  the  dependence  of  A^.  -  A^  q  on  6p  is  that  which  is  to  be  expected  from  the  fact  that  for  p  =  0, 
the  condition  of  equilibrium  is  that  A  be  a  minimum. 

There  will  be  a  first-order  effect  on  the  volume  of  the  system.  Let  be  the  volume  the  system 
would  have  at  zero  pressure  if  the  classical  statistics  held,  and  let  V  be  its  actual  volume.  A  classical 
system  at  volume  V  would  have  the  pressure  6p.  To  go  to  zero  pressure  its  volume  must  change  by 
“(0V/Bp)fTi^P  6p.  We  therefore  have 

V  -  Vg  =  (ev/ep)T^c^p  =  -RT(eT)/0v)rp  (ev/3p)T,c  (7) 

There  will  also  be  a  first-order  effect  on  the  energy.  We  have,  in  general,  for  the  sum  of  the 
thermal  energy  and  zero -point  energy 

Et  +  Eg  =  RT^  d  In  Q/dT  (8) 


*In  general,  we  use  a  subscript  c  to  designate  the  classical  value  of  any  quantity. 
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At  a  given  volume  V  therefore 

Et  +  Ez  -  Et,c  =  RT^  [(91n  Q/dT)y  -  (ain  Qc/aT)yj  =  RT^{dTi/&T)y  (9) 

where  ^  is  the  classical  thermal  energy  at  the  volume  V,  As  noted  above,  quantization  causes  no 
change  in  E  ,  since  any  effect  is  taken  care  of  by  17.  However,  at  volume  V  the  pressure  of  a  classical 
system  would,  as  we  have  seen,  be  6p.  Going  back  to  zero  pressure,  the  energy  would  change  by 
-(dE/aV)rj,  ^  (6V/ap)rp  c  6p.  Further,  at  zero  pressure  we  may  write  (0E/5V).ji  =  T(0S/5V),p  ^  = 

T(0p/0T)y^’^,.  Therefore,  for  fixed  zero  pressure  we  may  write  ’ 

E  -  Ec  =  RT2(aTj/eT)y  -  KT^(dT]/dV)^  (dp/dT)y^^  (dv/8p)^^^ 

=  RT2j^(072/0T)y  +  (aij/av)^  (0V/aT)p^g  j  (10) 

Provided  the  differential  coefficients  are  constant,  both  equations  7  and  10  are  exact. 

The  quantum  contribution  to  the  specific  heat,  ^Cp,  can  be  obtained  by  differentiating  equation  10. 
We  note  that  (d/dT)^  -  (d/dT)y  +  (0V/dT)p  (d/0V),p.  This  expression  is  to  be  applied  to  the  true  quan¬ 
tized  system  rather  than  a  hypothetical  classical  one,  so  the  (SV/^T)p,  which  appears  here  (without  sub¬ 
script  c)  is  the  actual,  rather  than  a  classical  value.  However,  the  operator,  (^/dT)p,  or  (d/dT)y  and 
{d/dV)rj^,  applied  to  a  classical  expression  gives  a  classical  expression.  We  therefore  obtain 

6Cp  =  2RT[(0??/aT)y  +  (ejj/ev)^  (ev/aT)p^j.  ] 

+  RT2|(d»]2/gT2)^  +  {d^ri/dTdV)  [id^'/8T)^  +  (0V/eT)p^p  ]  (11) 

+  {d\/dv^)y  (ev/eT)p  (av/aT)p^^  +  (STj/ev)^  (s^v/dT2)p^^j 

In  the  case  of  liquid  neon,  (0V/9T)p  and  (0V/9T)p^c  will  be  nearly  enough  alike  so  it  will  not  be  neces¬ 
sary  to  distinguish  between  them. 

2.  THE  EVALUATION  OF  THE  QUANTUM  CORRECTION 

The  evaluation  of  the  quantum  correction  for  the  solid  is  very  easy  if  we  assume  that  it  obeys  the 
Debye  specific  heat  law.  In  this  case  we  can  write  the  expression  for  In  Q  in  the  form^ 

In  Q  =  ln[(T/0)2e]  -  (3/4O)(0/T)2  +  (3/672 O)(0/T)4-  .  .  .  (12) 

where  B  is  the  Debye  characteristic  temperature.  If  we  break  off  equation  12  at  the  second  term  and 
compare  with  equation  1,  we  see  that  this  gives 

6  =  -(3/40)e2/3  Qg2/3=  -0.1464  Qc2/3  (13) 

In  many  respects  the  liquid  resembles  the  solid,  and  we  shall  assume  that  equation  13  holds  for 
the  liquid  as  well  as  the  solid.  It  is  true  that  there  is  a  translational,  as  ^ell  as  a  vibrational,  contri¬ 
bution  to  the  partition  function  of  the  liquid,  and  there  may  be  some  complications  arising  from  the  pos¬ 
sibility  of  a  change  in  coordination  number  in  the  liquid,  as  indicated  in  the  discussion  following  equation 

3.  However,  the  quantum  correction  will,  in  general,  be  expected  to  decrease  as  the  partition  function 
increases,  and  in  view  of  the  reasonably  close  similarity  between  liquid  and  solid,  it  would  appear  that 
equation  13  cannot  be  too  greatly  incorrect  for  the  liquid. 

We  may  now  refer  to  equation  2  to  get  V  in  terms  of  thermodynamic  quantities.  Since  A^,  -  Ep  = 

^c  "^^c  “  ^p  ”  ^t,c  “  "^^C)  where  E^^^.  is  the  classical  thermal  energy,  we  may  write  from  equation  13 

V  =  -0.146  e^®t,c/3RTg_2Sc/3R 

Since  (9Et^(./(9T)y  =  and  (5S(./0T)y  =  we  have 

(0?7/0T)y  =  -(2/3)7}Et^c/RT2 


(15) 
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We  alvSO  have 


{8v/dv)r^  =  (2/3)t)  |[OEt,c/av)T/RT]  -  [(as^/ev)x/R  ]  | 

CO 

rH 

We  note  that  T(aS^/av)q'  =  {dE^/dV)rjy  +  p.  Substituting  this  into  equation  16,  and  remembering  that  E^  = 

E^  ^  +  Ep  (the  zero-point  energy  does  not  occur  in  the  classical  case)  we  get 

(arj/o7)x  =-(2/3)rj  [^(dEp/dV)/RT  j  -  (2/3)  n  [p/Rt] 

(17) 

We  will,  in  general,  apply  this  with  p  =  0. 

To  apply  equation  11  we  also  need  the  second  derivatives.  In  the  final  expressions  for 
derivatives  we  set  p  =  0. 

the  second 

((P‘ri/eT\^  (4/9)77  [Et^g^/R^T^)  -  (2/3)77 

,, 

+  (4/3)77(Et^c/RT3j 

(18) 

d277/8VdT  =  (4/9)77  [E|.^g(dEp/dV)/R2T3] 

-(2/3)77  [OE.^^/0\7).p/RT2] 

(19) 

(e^77/av^).p  =  (4/9)77  [  (dEp/dV)2/R2T2  j 

-(2/3)77  [  (d2Ep/dv2)/RT  ] 

(20) 

-(2/3)n  [(ep/ev)^/RT] 

equation  19  can  be  somewhat  transformed.  Using  the  relation  at  zero  pressure, 

(eEc/0V)T  =  (SEt,c/8V5T  +  (dEp/dV)  =  T{8Sc/8V)j  =  T(0p/aT)v^c 

we  have 

8^r)/8V8T  =  (4/9)  rj  [Et^c(dEp/dV)/R2T3  ]  +  (2/3)n  [(dEp/dV)/RT2]  -  (2/3),,  [(0p/0T)y/RT]  (21) 

3.  APPLICATION  TO  SOLID  NEON 

Application  of  the  above  results  to  a  Debye  solid  should  yield  nothing  essentially  different  from  ap¬ 
plication  of  the  Gruneisen  equation,^ 

p  =r(Et  +  E2)-  (dEp/dV),  (22) 

where 

r=-0~ld0/dV  (23) 

However,  by  way  of  illustration,  it  will  be  of  some  interest  to  apply  equations  7  and  11  to  solid  neon.  In 
this  case,  we  shall  not  use  equations  15  or  17  but  will  go  directly  to  the  relation  which  follows  from 
equation  12. 

Tj=  -(3/4O)(0/T)2  ^  (3/6720) (0/T)4  (24) 

From  this  we  obtain 

(8v/8V)j.  =r[(6/4O)(0/T)2  -  (12/672O)(0/T)^  ]  (25) 

=  -2rrj[l  -  0.0060  (0/T)2  ] 

We  may  also  v/rite 

(av/ap)rp^^  =  ~(dV/dT)rj.^^/{^p/OT)y^^  (26) 

(c»p/aT)y  ^  is  readily  obtained  by  differentiating  equation  22,  taken  for  the  classical  case.  Since  F and 
dEp/dV  are  independent  of  T,  this  gives 

(0p/ST)v^c  =  Cj.,c  SRO 


(27) 
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(0V/dT)  may  be  written  in  the  form 
P>c 

(0V/0T)p^^  =  (dine /din  T)p^^  (28) 

Using  equations  25  to  28  in  equation  7  we  obtain 

r(V-Vc)  =  (2/3)(81n0/einT)p^p?7[l  -O.OO6O(0/T)2j  (29) 

If  we  have  a  series  of  substances  which  obey  the  law  of  corresponding  states,  the  value  of  (01n  0/ din  T)p^c 
will  be  the  same  for  all  of  them.®  In  the  case  of  the  rare  gases,  xenon  may  be  considered  as  practically 
classical  at  its  melting  point.  The  melting  points  are  very  near  to  corresponding  states;  hence  we  may 
take  the  value  of  (ain  S/dln  T)^  given  by  Murphy  and  Rice  for  xenon  at  its  melting  point*  as  the 
(a  In  @ /a In  T)p  c  for  neon  at  its  melting  point. 

If  we  set^  c  equal  to  the  value  of  0  at  V^,  then  it  is  seen  that  the  left-hand  side  of  equation  29  will 
be  (0-  ®c)/^av>  where  0^^  may  be  taken  as  the  mean  of  0  and  0^,  0  at  the  melting  point  of  Ne  (or 
rather  the  temperature  of  the  corresponding  point,  based  on  the  melting  point  of  xenon,  which  differs  by 
only  a  few  tenths  of  a  degree)  may  be  found  by  a  slight  extrapolation  in  Table  V  of  Murphy  and  Rice, 
while  0Q  may  be  estimated  from  the  results  with  xenon.  We  find 

(0-0  )/0  ,  =  (59.2  -  65.8)/62.5  =  -0.106  (30) 

c  a\ 

Evaluating  the  right-hand  side  of  equation  29,  using  0  =  59.2  and  T  =  24.9  to  calculate  r?,  we  obtain 

-(2/3)(ain0/01n  T)p^c^  [l  -  O.OO6O(0/T)2  ]  ^  .^o.667  x  0.385  X  0.410  X  0.965  =  -0.101  (31) 

This  is  an  excellent  check,  especially  since  (0V/0p)x,c  will  actually  vary  with  volume,  but  is  evaluated 
at  Vc. 

The  value  of  F  for  argon  has  been  taken  as  0.18  moles/cc.  For  neon  it  can  be  calculated  by  multi¬ 
plying  by  the  inverse  ratio  of  the  critical  volumes®  according  to  the  law  of  corresponding  states;  it 
therefore  is  0.325  moles/cc.  This  gives  V  -  Vp  =  0.33  cc/mole.  The  actual  volume  of  solid  neon  given 
by  Clusius'^  is  13.98  cc/mole.  If  we  multiply  the  volume  of  solid  krypton®  at  its  melting  point  by  the 
ratio  of  the  critical  volumes  we  get  13.42  cc/mole  for  V^.  Similarly,  from  argon®  we  get  13.59  cc/mole 
for  Vp.  The  calculated  value  of  V  -  seems  to  agree  about  as  well  as  one  might  expect. 

We  have  also  evaluated  equation  11  for  the  solid.  This  requires  some  further  differentiation,  in¬ 
cluding  differentiation  of  equation  28,  which  brings  in  the  quantity  (a^ln  0/dln  T  0T)p^p,  which  has  been 
estimated  as  -0.00405  X  (289.8/44.8)  from  Table  H,  Murphy  and  Rice  (the  ratio  in  parentheses  is  the 
ratio  of  critical  temperatures  of  xenon  and  neon,  taken  from  Guggenheim).  We  also  have  to  evaluate 
(0V/dT)p  (nonclassical)  which  is  done  with  the  aid  of  Table  V  of  Murphy  and  Rice.  We  have  obtained  for 
«Cp  at  the  melting  point  2.0  cal/mole/deg.  Comparison  of  Tables  II  and  V  of  Murphy  and  Rice  shows 
that  this  is  about  right. 

4.  APPLICATION  TO  LIQUID  NEON 

The  checks  obtained  in  Section  3  merely  confirm  the  consistency  of  the  present  calculation  with  that 
of  Murphy  and  Rice.  The  comparison  of  experimental  and  calculated  V-  may  have  some  significance, 
but  again  the  calculation  may  be  handled  equally  well  by  either  method.  For  the  case  of  the  liquid,  how¬ 
ever,  only  the  equations  of  the  present  paper  are  available. 

To  calculate  V  -  Vp,  or  better  (V  -  Vp)/Vp,  we  clearly  need,  from  equations  7,  17,  and  26,  values 
of  Vp,  dEp/dV,  (aV/aT)p^P  and  (ap/0T)y^p.  To  get  6Cp  we  need  also  E^^p,  Cj,^p,  (d^V/dT^)^^^  and  d^E/ 
dV^.  These  quantities  are  not  known  for  neon;  but  they  are  known  or  may  be  readily  calculated  for  argon 
from  related  quantities.^®  Since  the  quantities  appearing  in  the  expressions  for  (V  -  Vp)/Vp  and  ^ Cp  are, 
except  for  17,  in  a  combination  which  is  independent  of  the  substance  if  the  law  of  corresponding  states 
holds,  the  values  for  argon  should  be  sufficiently  good.  The  values  for  the  second  derivatives  show  a 
rather  sudden  change  at  a  molar  volume  equal  to  the  molar  volume  at  the  normal  melting  point  of  argon. 
However,  the  volumes  we  are  interested  in  are  equal  to  or  greater  than  the  normal  melting  volume,  so 
no  ambiguity  occurs  here.  The  second  derivatives  can  be  obtained  by  differences  from  the  first^ deriva¬ 
tives;  the  values  are  necessarily  rough,  but  they  fortunately  do  not  contribute  greatly  to  ^Cp,  d  Ep/dV^ 


*We  will  use  the  calculations  based  on  Murphy’s  and  Rice’s  Curve  II.  These  calculations  do  not  fit 
the  experimental  data  on  neon  and  xenon  perfectly,  but  they  are  reasonably  good,  and  they  do  exactly 
fit  the  law  of  corresponding  states  as  extended  by  Murphy  and  Rice. 
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is  sufficiently  small  and  the  terms  containing  it  are  sufficiently  unimportant  that  they  can  be  neglected. 

We  have  yet  to  consider  r]  for  neon.  This  is  to  be  obtained  from  equation  13.  can  be  obtained 
by  comparison  with  argon.  We  have  previously  shown  that  can  be  expressed  approximately  in  the 
form^° 

Qc  =  Q§  =[(2e)^/3(T/e)  +  (27rrakT/h2)l/22-l/6e(a  -  ao)]^  (32) 

On  going  from  one  substance  to  another  which  obeys  the  law  of  corresponding  states,  the  two  terms  in 
the  bracket  transform  in  the  sam,e  v/ay.  (For  the  transformation  of  &  see  Murphy  and  Rice.®  To  con¬ 
vert  e  from  argon  to  neon,  we  use  (argon)  =  1.530  and  f^  (neon)  =  1.405.)  Qq  has  been  calculated 
for  argon  in  connection  with  the  calculations  of  the  article  cited,  v/hence  Qq  for  neon  and  17  for  neon 
can  be  obtained  as  indicated.  The  various  thermodynamic  quantities  required  and  r}  are  listed  in  Table 
1, 


Table  1.  Thermodynamic  Quantities  for  Argon;  7]  for  Argon  and  Neon. 


V  =  28.02  cc/mole 
T(fusion)  =  83.78  °K 


dEp/dV  =  50.5  cal/cc 
Et/3RT  =  0.7961 
Cj,/3R  -  0.856 

(aV/0T)p  =  0.1178  cc/mole/deg 
{dp/d'T)y  =  0.556  cal/cc/deg 
(d^V/dT^)p  =  0,00165  cc/mole/deg^ 
QQ(argon)  =  3.681 
?7  (argon)  =  0,01080 
V  (neon)  =  0.1030 


This  table  gives  the  values  of  the  various  quantities  at  the  melting  temperature  and  melting  volume 
of  argon,  or,  in  the  case  of  r?,  at  the  corresponding  point  for  neon.  Neon  at  its  melting  point  (24. 55°^), 
or  at  its  corresponding  temperature  (24, 9® K,  which  is  the  temperature  we  have  actually  used),  is 
slightly  expanded  beyond  its  corresponding  volume.  The  errors  which  will  be  caused  by  this  discrepancy 
are  not  sufficient  to  make  it  worthwhile  to  perform  further  calculations.  It  is  easy  to  find  from  Refer¬ 
ence  10  how  the  various  quantities  involved  change  with  volume  at  constant  pressure.  This  should  give 
the  order  of  magnitude  of  the  errors  involved,  and  it  is  found  to  be  of  the  order  of  20%  in  (V--  Vc)/Vc 
and  20%  in  6Cp.  However,  these  are  probably  exaggerations,  since  the  most-affected  quantity  is  tj,  and 
it  will  be  affected  considerably  less  at  constant  temperature  than  at  constant  pressure. 

Use  of  Table  1  enables  us  to  estimate  (V  —  V^)/V^  as  0.026.  Figure  2  of  Guggenheim’s  article®  en¬ 
ables  one  to  estimate  (V  “  V^)/V^  by  comparing  the  ratio  of  density  to  critical  density  for  neon  and 
other  gases.  It  appears  that  neon  is  about  4  to  5%  expanded  over  what  one  might  expect  from  the  law  of 
corresponding  states.  However,  just  the  points  near  the  melting  points  near  the  melting  point  for  neon 
are  a  little  out  of  line,  so  4%  expansion  is  probably  the  better  estimate.  Pitzer  states  that  liquid  neon 
has  a  density  of  about  3%  less  than  that  of  the  "'perfect  liquid. So  it  appears  that  the  theoretical 
estimate  of  (V  -  V(,)/V(.  cannot  be  far  out  of  line. 

6Cp  is  calculated  as  1.53  cal/m,ole.  The  specific  heat  of  liquid  neon  is  given'^  as  8.64  near  the  melt¬ 
ing  point,  that  of  argon"^  as  10.05  or^^  10,5,  that  of  krypton^®  as  10.6,  and  that  of  xenon^^  as  10.7. 

Therefore,  it  is  obvious  that  there  is  reasonably  good  agreement  between  the  calculated  and  exper¬ 
imental  values. 
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5.  TWO-PHASE  SYSTEMS 

We  will  now  consider  the  effect  of  the  quantum  correction  on  the  equilibrium  between  the  solid  and 
the  liquid.  We  shall  first  consider  the  change  in  the  temperature  of  melting  at  zero  pressure.  At  zero 
pressure  the  Gibbs  free  energy  F  is  equal  to  the  Helmholtz  free  energy  A.  We  can  consequently  use 
equation  3,  substituting  F  for  A  and  F^  for  A^..  F^,  really  corresponds  to  A^,  q,  since  we  should  compare 
classical  and  quantum  F’s  at  the  same  pressure.  The  difference  between  A^,  and  A^^o  is,  as  we  have 
shown  following  equation  6,  a  second-order  quantity.  It  is  not  entirely  negligible  for  the  solid,  but  is 
not  important  and  may  be  neglected  for  our  purposes. 

The  condition  for  equilibrium  between  two  phases  may  be  written 

AF  =  -  RT  An  =  0  (33) 

If  there  were  no  quantization,  equilibrium  would  occur  at  a  temperature  T^  with  Ap^  =  0,  while  equation 
33  holds  at  a  temperature  T  where  A  F^.  =  0  +  (dAF^,/ST)p  (T  -  T^.)  =  -^S^{T  ~  Substituting  this  into 
equation  33  we  obtain 

T  -  T^  =  -RT  An/  A^c  (34) 

Using  the  values  of  n  for  solid  and  liquid  given  in  Sections  3  and  4,  and  taking  for  AS^,  the  entropy  of 
melting  of  argon,  3.35,  we  obtain  T  “  T^,  =  -4.4*".  The  observed  value  of  T  -  T^  is  about  24.55  -  24.9  = 
-0.35®,  The  discrepancy  is  very  large  and  quite  surprising,  considering  the  the  apparently  good  results 
obtained  with  solid  and  liquid  separately. 

We  can  analyze  the  situation  further  by  considering  AE  ~  AE^  and  AS  —  AS^.  for  the  fusion  process. 
Since  the  melting  temperature  is  so  nearly  a  corresponding  temperature  with  respect  to  the  melting 
temperature  of  xenon,  say,  we  can  neglect  any  effect  of  change  in  temperature  on  these  quantities.  This 
means  that  we  can  calculate  aE  ~A  E^  by  applying  equation  10  directly  to  both  solid  and  liquid  at  24.9® K. 
For  the  solid,  we  apply  the  equations  of  Section  3.  For  the  liquid,  we  apply  the  equations  of  Section  2  and 
the  data  of  Section  4,  and  convert  from  argon  to  neon.  We  obtain  E  -  E^  =  -  54.3  cal/mole  for  solid  argon 
and  E  -  Eg  =  18.3  for  liquid  argon.  Thus, 

AE  “AEg  =  -36.0  cal/mole,  calc. 

From  the  heat  of  krypton,  we  get  for  the  corresponding  heat  of  fusion  of  neon  83.2  cal/mole,  while  the 
actual  heat  of  fusion  is  about  80.1  cal/mole.  Thus, 

AE  -  AEg  =  -3.1  cal/mole,  obs. 

A  similar  discrepancy  occurs  for  the  entropy.  We  can  write 

AS  -  ASg  =  (AE  -  AEg  -  AA  +  AAg)/T 

=  [(AE -AEc)/T  ]+ RAtj 

=  -1.45  -}-  0.61  =  -0.84  cal/mole/deg 

The  actual  value  of  AS  “ASg  is  around  -0.1  cal/mole/deg. 

The  results,  with  AE  “AEg  and  AS  “  AS^,  appear  to  indicate  that  the  quantum  correction  for  the 
liquid  is  nearly  as  great  as  for  the  solid.  This  would  seem  very  strange  in  itself,  but  it  is  still  more 
puzzling  in  the  light  of  the  results  of  Section  5.  The  calculation  for  6Cp  may  indeed  indicate  that  equa¬ 
tion  13  gives  something  of  an  underestimate  of  t];  if  it  were  increased  by  perhaps  30%,  a  better  value 
of  Cp  would  result.  On  the  other  hand,  it  must  also  be  noted  that  the  derivatives  of  tj  in  equation  11  in¬ 
volve  the  quantity  dEp/dV,  which  is  somewhat  indirectly  calculated  from  the  data.  It  does  not  appear 
likely,  however,  that  this  can  cause  appreciable  error.  We  have 

dEp/dV  =  (aE/aV)^  “  idE^/dYjry, 

(dEt/8V)'p  is  certainly  negative;  therefore  dEp/dV  should  be  slightly  larger  than  (aE/aV).^  which  at  zero 
pressure  is  equal  to  T(8p/dT)y.  Now  the  value  of  dEp/dV  used  is  only  slightly  larger  than  TOp/0T)y; 
it  seems  likely  that  dEp/dV  should  be  increased  rather  than  decreased,  and  this  would  actually  increase 
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I^CpI  ,  so  that  |t|  I  would  have  to  be  increased  less.  The  use  of  values  at  the  normal  melting  temperatures 
and  volumes  gives  a  slight  overestimate  for  |6Cp|,  which  might  mean  that  |t]|  should  be  increased  by  a 
little  more  than  30%,  but  this  effect  is  not  very  great. 

Another  quantity  that  occurs  in  the  derivatives  of  tj  that  is  not  entirely  directly  determined  experi¬ 
mentally  is  We  have  also  used  the  calculated,  rather  than  the  experimentally  determined,  value 
of  Cp  However,  it  does  not  appear  that  any  possible  change  in  these  quantities  could  clear  up  the 
discrepancy. 

On  the  other  hand,  it  appears  likely  that  the  discrepancy  is  too  great  to  be  accounted  for  on  the 
basis  of  deviations  from  the  law  of  corresponding  states. 

Therefore,  we  are  forced  to  leave  the  question  of  the  quantum  correction  for  neon  in  a  somewhat 
unsatisfactory  state. 
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